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We report on our lattice calculations of the nucleon's generalized parton distributions (GPDs), 
concentrating on their first moments for the case of Nf = 2. Due to recent progress on the nu- 
merical side we are able to present results for the generalized form factors at pion masses as 
low as 260 MeV. We perform a fit to one-loop covariant baryon chiral perturbation theory with 
encouraging results. 
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1. Introduction 

The study of the internal structure of hadrons still presents an exciting challenge. Among the 
different types of studies, the computation of Generalized Parton Distributions (GPDs) is especially 
challenging, but also attractive, because of their potential for hadron physics. 

GPDs were introduced in the late 90s. For a given hadron, they provide detailed information 
on the partonic structure with respect to spatial, momentum and spin degrees of freedom. GPDs 
combine the information of the traditional form factors and parton distribution functions (contain- 
ing them as limiting cases) into a single set of functions and hence contain information also on 
the correlation between the momentum, spin and spatial degrees of freedom. For the nucleon, one 
hopes GPDs will provide one day a three-dimensional spatial picture, a better understanding of its 
spin structure and a value for the quark orbital angular momentum [1]. 

Beside the renormalization scale 1 , GPDs depend on three kinematic variables: the longitudinal 
parton momentum fraction x, the skewness parameter t, and the virtuality t. The quark structure of 
a nucleon, for example, is governed by eight GPDs. Among these, the unpolarized GPDs H and E 
parametrize the off-diagonal matrix element 

(N(P')\^(x)\N(P))=U(P')^H(x^,t) + ia ^E(x,^t)\u(P) +highertwist. (1.1) 

Here P and P' denote the incoming and outgoing nucleon momenta (and so A = P' — P, P = 

(P'+P)/2, t = A 2 and t, = -n ■ A/2) and 6§ (x) is the light-cone bilocal operator 



M ,ikx -,( * \ ,Ai ( y i) -ig^ ,1 danA(an) A 



e'^ql- ~ n \t&e- l8i W aanA[an >q i-nj , (1.2) 

which often arises in hard scattering processes (see, e.g., [3]). n (with P -n = 1) denotes a light cone 
vector in Eq.(1.2) and 3? the correct path-ordering of the gluon fields A. The polarized nucleon 
GPDs, H and E, are defined in a similar manner, replacing y^ in Eq.(1.2) by 7^/5. 

2. GPDs and the lattice 

GPDs can be accessed experimentally, for instance, via deeply virtual Compton scattering. 
The analysis, however, is rather demanding and requires also a partial modeling of the combined 
x-, <§- and ^-dependence. Cross-checks to other methods are thus inevitable. 

A promising method is given by lattice QCD computations. Although a direct determination 
of GPDs on the lattice is not possible, their (Mellin) moments 

I dxx n - l H(x,£,,t), J dxx n ' l E{x^,t), ... (2.1) 

are accessible. For a nucleon, for example, these moments can be calculated via matrix elements 
(N(P')\0\N(P)} of local operators O. For H and E these operators read 

0$ V[V "- 1 (z) = q{z) r^iD Vl • • • iD v - i} q{z) - traces (2.2) 



'For simplicity, we drop the explicit reference to the renormalization scale /i in what follows. It is always implicitly 
understood. Our lattice data below has been nonperturbatively renormalized [2] and is for the MS scheme at fl = 2GeV. 
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— — 
where q refers to a quark field, D = D — D to the covariant derivative and { • • • } to a symmetrization 

of the Lorentz indices. For a definition and further details on the operators needed for the remaining 

nucleon GPDs, the reader may refer to [4]. 

Admittedly, the computation of such matrix elements is quite demanding already for n>2, 
and we are not yet in the stage to provide precision results close to the physical point. Nonetheless, 
such calculations have become more and more feasible in recent years, and hence have attracted 
interest from within the lattice community [5, 6, 7, 8, 9, 10, 11]. 

In what follows, we will restrict ourselves to the two nucleon GPDs H and E. Their moments 
are polynomials in E, , 



dxx' 



n-\ 






[(»-l)/2] 
I (2§) 

k=0 



2k 



A n ,2k(t) 
B n ,2k(t) 



±S n , ev6a (2$) n C n {t). 



(2.3) 



The expansion coefficients A, B and C are real functions of the momentum transfer t (and the 
renormalization scale /i) and are known as the Generalized Form Factors (GFFs) of the nucleon. 
In this notation, for instance, A\$ and 5io correspond to the electromagnetic form factors [12], and 
^20> #20 an d ^2 parametrize the matrix elements of the energy-momentum tensor Oy V 



{N(P f 



K V \N(P)) 



U{P')\y^P v] A 20 (t) 



iA p oP&- 

IfflA 



P V} B 20 {t) 



A^A V > 



m N 



C 2 (t)}U(P). (2.4) 



Below we present results for A 2 q, B 2 q and C 2 . They can be extracted from ratios of two- and 
three -point correlation functions 



R(t,X,p',p) 



C3(t,T,p',p) 
C 2 (t,p') 



C 2 (z,p')C 2 (t,p')C 2 (t-x,p) 
C 2 (z,p)C 2 (t,p)C 2 (t-z,p') 



1/2 



(2.5) 



which are proportional to (N(P')\Oy V \N(P)) and constant in the limit < x < t < T/2 (T is the 
temporal lattice extension). C 2 (t,p) is the nucleon two-point function with a source at time and 
sink at time t, and C3 (t , X, p' ,p) is the three -point function with an operator insertion at time x. The 
latter we calculate employing the sequential source technique. 



3. Results 
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Our data for the GFFs is for gauge configura- 
tions thermalized using the standard Wilson gauge 
action and two flavors of clover-improved Wilson 
fermions. The gauge couplings are j3 = 5.25, 5.29 
and 5.40; and the K values are such that pion masses 
from 1 GeV down to 260 MeV are simulated, where 
we primarily work with the data in the mass range 
260 MeV < m n < 490 MeV. The scale is fixed 
through setting ro = 0.5 fm. This is about the value 
we obtain from chiral extrapolations of our nu- 
cleon mass data [13] for the same set of config- 
urations. The lattice sizes are 24 3 x 48, 32 3 x 64, 40 3 x 64 and 48 3 x 64. In particular the latter 
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Figure 1: A 2 o(t) for the isovector case for /3 = 
5.29 and K = 0.13632 and for different volumes. 
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Figure 2: The generalized form factors A20, B20 and C2 (from top to bottom) vs. momentum transfer —f, 
left for the isovector channel, right for the isosinglet channel. The data is for three lattice spacings and two 
groups of approximately equal pion masses. If applicable, solid (dashed-dotted) lines represent dipole fits 
to the data. Dashed lines at low t result from a simultaneous fit of the low-pion-mass data (full symbols) to 
covariant chiral perturbation theory (see text and also Fig. 3 for more details). 
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two provide us with a good signal-to-noise ratio. See, for example, Fig. 1, where data for A%Q{t) 
in the isovector channel is shown for the lattice sizes 24 3 x 48, 32 3 x 64 and 40 3 x 64 at /3 = 5.29 
and K = 0.13632 {m n = 287 MeV). The number of measurements is 2755, 3013 and 1478, respec- 
tively. The advantage of volume averaging is clearly evident as with about half the statistics, the 
data for the 40 3 x 64 lattice comes with much less statistical noise than that for the 24 3 x 48 lattice. 
Moreover, Fig. 1 indicates that finite size effects are small, at least at our level of precision. 
A selection of all of our GFF data is shown 
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as a function of the momentum transfer —t in 
Fig. 2. There, the panels from top to bottom dis- 
play the respective data for A2o(t), #20 (0 an d 
C2(t). Left panels are for the isovector case, right 
panels for the isoscalar (without disconnected con- 
tributions). For simplicity, only data for five en- 
sembles is shown, which fall into two groups of 
approximately equal pion mass. For the larger 
pion mass (430-490 MeV), we have results for 
three lattice spacings (a = 0.06, 0.07 and 0.08 fm), 
for the smaller one (260-287 MeV) we can show 
data for two sets (a = 0.06, 0.07 fm). For this 
(admittedly small) range of lattice spacings we 
observe, however, no systematic dependence on 
a. Apparently, there is a slight vertical shift in 
the data for A2o(t < 0) [and in the opposite direc- 
tion for #20 (/ < 0)] f° r the lighter sets, but we do 
not see these shifts for the heavier sets (including 
those not shown), at least with the available pre- 
cision. It will be interesting to see how well our forthcoming results at j8 
a = 0.084 fm, 260 MeV pion mass) fit to these findings. 

Similarly, we observe a trend for A2o(t) if the pion mass is changed: The low-? dependence 
of A 2 o(t) gains slope if m K is reduced from 430-490 MeV to 260-287 MeV. This effect, however, 
is small and we see no such effect in the data from 500 MeV to 1 GeV pion mass. It thus remains 
to be seen if this effect at lower m % stays or disappears with higher statistics. As above, a further 
check should become possible as soon as our results at j8 = 5.25 and K = 0.13620 are available. 

We can confirm though the (notorious) weak m n dependence of A u ^ d at t = 0, i.e., of (x) u ~ d 
(see upper left panel in Fig. 2). From phenomenology one expects (x) u ~ d m 0.16 at the physical 
point. So far, however, all available (world) lattice data for (x) u ~ d for pion masses above 200 MeV 
gives values for (x)"~ d well above 0.16, and moreover, almost no signal for a downwards trend 
towards the physical point is seen (see, e.g., [4] and references therein). From baryon chiral pertur- 
bation theory (BChPT), for example, such a trend is expected, but it has not yet been demonstrated 
(convincingly) on the lattice. 

It is however interesting that our data for \t\ < 0.4GeV 2 indicates an almost linear ^-dependence 
for A l 2Q d (t) and a flattening of the slope for B^q (t). This would be consistent with expectations 
from covariant BChPT at leading-one-loop order [14]. 



Figure 3: A^ a (top) and B 1 ^" (bottom) vs. -t at 
the pion masses 261 and 288 MeV. Dashed lines re- 
sult from a simultaneous fit of the data (including 
that for C^ d ) to expectations from BChPT [14]. 
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In Fig. 2, and in particular in Fig. 3, we show a first attempt of fitting our data to the BChPT 
expressions for A20, 620 an d C2 as worked out in [14] (see the dashed lines at lower t). Note 
that such a fit has to be a combined fit to the data for all three GFFs simultaneously, because 
the parameter a.20 enters all of them. The dashed lines in Figs. 2 and 3 refer to such a fit which 
incorporates only the lighter data sets (full symbols) and points for \t\ < 0.44 GeV 2 . Five parameters 
(«20> &20 C20, Cg and cu) were left free, while the phenomenological value (Ax)jj _™ = 0.21 was 
used to constrain the coupling Aa^g. The m^-dependence of the nucleon mass (entering the BChPT 
expression for £20 and C%) and the parameter Mq were taken from our nucleon mass fits [13]. 

It turns out that for the isovector case the fit quality is quite good: A reduced # 2 -value of about 
one is reached and the low-? dependence of A20 and B20 is roughly reproduced; actually, also for 
C%, as only the data point at the smallest \t\ falls somewhat below the fitting curve. For the isoscalar 
case, the same fit works less satisfactorily. The reason might be that BChPT does not work at the 
pion masses under consideration, or that disconnected contributions are not negligible. The latter 
are certainly worth to be calculated. 

Even though disconnected contributions are still 
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J" : = 5.29 H 
4 P = 5A0 I- 



missing, it is interesting to look at the total quark °- 4 " " j"-. = 5.29 ^ 

angular momentum 0.3 - tT £ • | = j^-j j: 

0.2 - 
(3.1) 



ju,d _ 



2 



A^(0)+<(0) 



0.1 
0.0 



. ~--U 

to check if it is in the ballpark of expected values. 

In Fig. 4 we show this data for our lighter sets, that J 

is, for pion masses 261 and 287 Me V. Note that °' 2 [ 

for Ajq d we have data directly at t = 0, but not ' m [GcV , 

for #2^(0). However, looking at Fig. 2 one easily 

sees that the main contribution to J u+d comes from Figure 4: The total quark angular momentum J 

A l £ d (0) and the ^-dependence of B'£ d is compa- and s P in s '' vs - P ion mass - 

rably weak. We therefore approximate B^q (0) by 

our data for the smallest \t\. This should be perfectly fine for our purposes, given all the other 

uncertainties and the lack of disconnected contributions. Note that in Fig. 4 we have also included 

data for the quark spin j 

s q = -J_ i dxH(x^,0) = -A\ (t = 0), (3.2) 



which we obtain from data for the axial nucleon GFFs A" and A "j, . 

If we compare our data in Fig. 4 with that of the LHPC collaboration [15], we find good 
agreement (albeit their data is for Nf = 2 + 1). We also see the same ordering for the total and 
orbital {L q = J q — s q ) angular momentum and the quark spin: 

\J d \ < \J% \J d \ < |/|, \L u+d \ < \L u \,\L d \ . 

It will be interesting to see how this figure changes when data at smaller pion masses becomes 
available and/or disconnected contributions are included. 



^TTnfnctunotoK/ tliavo ic r\r\ Hoto £r\r il"^" 



Unfortunately, there is no data for AjJ at f5 = 5.4 for these small pion masses, but it will become available soon. 
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4. Conclusions 

We have presented an update on our efforts to calculate the generalized form factors for the 
nucleon. We have restricted ourselves here to the case of Nf = 2 and reported only on results for 
the GFFs of the energy-momentum tensor (n = 2). Due to recent progress on the numerical side we 
are able to provide data for these GFFs for pion masses down to 260 MeV. In particular our lighter 
sets provide an improvement of the available data for these form factors: Large lattice volumes 
have allowed us to obtain a very good signal-to-noise ratio, and at low \t\ our data starts to fulfill 
expectations from one-loop BChPT. When comparing our data to that of the LHPC collaboration 
presented at this conference [15] we see a small vertical offset for the GFF data, but overall agree- 
ment for angular momentum and spin. It remains to be seen if this offset is due to the different 
renormalization procedures of the lattice operators or due to the different Nf. 

The numerical calculations have been performed on the APEmille, apeNEXT systems and PAX cluster 
at NIC / DESY (Zeuthen, Germany), the IBM BlueGene/L at EPCC (Edinburgh, UK), the BlueGene/P 
(JuGene) and the Nehalem Cluster (JuRoPa) at NIC (Mich, Germany), and the SGI Altix and ICE 8200 
systems at LRZ (Munich, Germany) and HLRN (Berlin-Hannover, Germany). We have made use of the 
Chroma software suite [16]. The BlueGene codes were optimised with Bagel [17]. This work has been 
supported in part by the DFG (SFB/TR 55, Hadron Physics from Lattice QCD) and the EU under grants 
238353 (ITN STRONGnet) and 227431 (HadronPhysics2). A.St acknowledges support by the European 
Reintegration Grant (FP7-PEOPLE-2009-RG, No.256594). JZ is supported by the University of Adelaide 
and the Australian Research Council through a Future Fellowship (FT 100 100005). 
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